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Two of the most important and pressing questions in cosmology and particle physics 
are: (1) What is the nature of cold dark matter? and (2) Will near-future experiments 
on neutrinoless double beta decay be able to ascertain that the neutrino is a Majorana 
particle, i.e. its own antiparticle? We show that these two seemingly unrelated issues 
are intimately connected if neutrinos acquire mass only because of their interactions 
with dark matter. 



The existence of cold dark matter in the Universe is now well accepted pQ. From the 
viewpoint of particle physics, it should consist of a weakly interacting yet-to-be-discovered 
neutral stable fermion or boson. A prime candidate is the lightest supersymmetric particle 
(LSP) in the minimal supersymmetric standard model (MSSM). More generally, we need 
only an exactly conserved Zi symmetry EJ and some new particles which are odd under 
it, keeping all known particles even. In the MSSM, this Z 2 symmetry is called R parity, and 
the new particles are the squarks, sleptons, gauginos, and higgsinos. 

Consider now the interactions of the neutrino with particles in this new class. To realize 
the well-known dimension-five operator for Majorana neutrino mass j3j, 

£eff = ^f K0° - la<f> + )M° ~ h<P + ) + H.C., (1) 

where (v a ,l a ) and $ = (0 + ,0°) are the usual lepton and Higgs doublets of the standard 
model (SM), it is clear that the new particles must form a loop with four external legs given 
by VaVpcfPcfP. There are generically three such one-loop diagrams |3]. In the MSSM, this 
does not happen because this operator also requires lepton number to change by two units. 
However, if a neutral singlet superfield N is added, then Figure 1 is generated as a radiative 
contribution to the neutrino mass. 
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Figure 1: One- loop radiative neutrino mass in supersymmetry. 

We note that the particles in the loop all have odd R parity. Of course, in the context of 
supersymmetry, this also implies that v couples to iV directly through 0°, so that a Dirac 



mass already appears, and together with the heavy Majorana mass of N, the famous seesaw 
mechanism [6 j allows v to obtain a tree-level Majorana mass. That is why Figure 1 is always 
negligible in the MSSM with the addition of N, However, in a more general framework, 
the new particles in the loop may be the only source of neutrino mass |H] , and in that 
case there will be interesting phenomenological implications on lepton flavor transitions and 
neutrinoless double beta decay, as shown below. 

Basically, the argument goes as follows. In order for the new particles in the loop to be 
identified as cold dark matter with the correct value of their relic density in the Universe at 
present, their interactions with neutrinos and charged leptons must not be too weak. On the 
other hand, they are also responsible for the masses of neutrinos and their observed mixing 
in neutrino oscillations. This implies necessarily flavor changing transitions such as /x — > ej. 
In order to suppress the latter, the parameter space of neutrino masses is limited, thereby 
enforcing a lower bound on neutrinoless double beta decay. 

Of the three generic one-loop diagrams giving rise to a radiative neutrino mass, the 
simplest in terms of new particle content is given in Ref. jH]. The standard model is extended 
by adding three neutral singlet fermions N{ and a second scalar doublet (i] + ,f]°), which are 
odd under an exactly conserved Zi symmetry, keeping all SM particles even. In that case, 
the analog of Figure 1 is Figure 2, as depicted below. 
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Figure 2: One-loop radiative neutrino mass in the model of Ref. jH]. 



We note again that the particles in the loop are all odd, and that lepton number changes 
by two units as in Figure 1. The Z2 invariant Higgs potential is given by 

V = m]$ , $ + m2f| t ?]4- -Ai($ f $) 2 + -A 2 (r/ f ?7) 2 + A 3 ($ t $)(?7 t ?7) 

+ A 4 ($tr ? )(r ? t$) + ^A 5 [($t??) 2 + H.c], (2) 

with (0°) = v and (77°) = 0. Let us choose the bases where Ni and l a , l c a are diagonal, and 
consider the interactions of (v a , l a ) with JVj and (?7 + , 77 ) , i.e. 



'JV 



h ai (v a r)°-laV + )Ni + H.c. (3) 



Since (r/°) = is required to preserve the exact Z2 symmetry, there are no Dirac masses 
linking u a with Ni. In other words, even though iV* have heavy Majorana masses, the 
canonical seesaw mechanism is not operative. Further, the lightest among the new particles 
will be stable and becomes an excellent candidate for the cold dark matter of the Universe. 
We see thus that neutrinos acquire mass here only because of their interactions with dark 
matter. 

We note that Figure 2 depends on the existence of the A5 coupling of Eq. (J2J). If it were 
zero, we could assign the exactly conserved additive lepton number —1 to (n + ,i]°) and to 
Ni, in which case the neutrinos would stay massless. This means that it is natural for A5 
to be very small, which we will assume from here on. Without loss of generality, A5 may 
also be chosen to be real. We now define rf = {jin + irji)/y2 and obtain m\ — mj = 2X^v 2 , 
where mn(mi) is the mass of r)n(r)i). Using m\ = {m 2 R + mj)/2, the radiative neutrino mass 
matrix is then given by [S] 

(KA N ^U^f^ , , 

where 

I(x)='-i±r[-—\ I ^ • 
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Assuming that atmospheric neutrino mixing is maximal, the neutrino mixing matrix U 

which diagonalizes Ai u can be written as U = UP, where U is approximately given by 

/ C12 S12 s 13 e^ 5 \ 

U ~ -s 12 /V2H-c 12Sl3 e*VV2 c 12 /V2 + s 12 s 13 e iS /V2 -1/V2 , (6) 

V -su/y/2 - c 12 s 13 e i5 /V2 c 12 /^2 - s 12 s i3 e iS /y/2 1/^2 j 

/ e iai/2 q q\ 

p = o e iQ2 / 2 , (7) 

V l) 

and C12 = 008 6*12, S12 = sin#i 2 , s\ 3 = sin 6*i 3 , with tan 2 6*i 2 — 0.45, and si 3 ^ 0.2. 

We now assume the lightest iV to be dark matter. Call it N k . We need to calculate its 
relic density as a function of its interaction strengths h ak , its mass M k , and the masses of 
r) , rjR, and rji, which we take for simplicity to be all given by mo, with m > M k . Our goal 
is to obtain the observed dark-matter relic density of Qdh 2 — 0.12 PlllOj. 

The thermally averaged cross section for the annihilation of two N k s into two leptons 
is computed by expanding the corresponding relativistic cross section a in powers of their 
relative velocity and keeping only the first two terms. Using the result of Ref. [TT], and 
recognizing that lepton masses are very small, we have 

y y k (l-2r k + 2rl) 



(av) = a + b k v + 



a = 0, 



24ttM 2 



where 



r k = Mlliml + M 2 ), y\ = £ \h ak h 

a/3 



* |2 
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(9) 



Following Ref. [T2"] . the relic density of N k is then given by VL<ih 2 = Y^soMkh 2 / p c , where 
Yoo is the asymptotic value of the ratio n Nk / s, with Y^ 1 = 0.2Q4gJ M P \M k (3b k /x 2 ), s = 
2970/cm 3 is the entropy density at present, p c = 3H 2 /8ttG = 1.05 x 10 -5 /i 2 GeV/cm 3 is 
the critical density, h is the dimensionless Hubble parameter, Mp\ = 1.22 x 10 19 GeV is 
the Planck energy, and g* is the number of effectively massless degrees of freedom at the 



freeze-out temperature. Further, x/ is the ratio M k /T at the freeze-out temperature and is 
given by 

0.0764Mpi(66 fc /x/)c(2 + c)M k 



xj = In 



G?**/) 1 / 2 



1/2 



Using gj = 10 and c = 1/2 as in Ref. [12], we obtain 



M fr 



GeV 

h 



5.86 x 10~ 8 x/ /2 e x f 



Sl d h 2 
0.12 



(GeV)" 



2.44 x 10" 11 x) 



0.12 



fi> 2 



(10) 

(11) 
(12) 



where b k and y& are given in Eqs. (jSJ) and (0). Since h k jy\ is a function of M k and m , 
Eqs. (|TT|) and (JT2J) allow us to calculate Mfc and m in units of GeV for a given set of yl, Xf 
and Qdh 2 . 
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Figure 3: M k versus m /yk for y fc = 0.3,0.5,0.7, 1.0 (left to right) for VL d h 2 = 0.12, where y^ is defined in 
Eq. JSJ. 

In Figure El we plot M k versus m /y k for y k = 0.3,0.5,0.7, 1.0. As we can see from the 
figure, the dark matter constraint requires that M k increases as mo increases and for each 
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value of mo, it may be as large as rriQ. We also see that mo/yk cannot exceed 350 GeV or so 
in the perturbative regime yk & 1. This is a very strong constraint, because m /yk sets the 
scale also for lepton flavor transitions such as ft — > ej and the experimental upper bound of 
its branching fraction cannot be satisfied, unless some cancellation mechanism is at work. 



The branching fraction of \x — > ej is given in this model by ^3] 

i 



£(//-> ery) 



3a 



Q4n(G F m 2 ) 



2\2 



c 4 



30 GeV 

, m o/C . 



where 



C 2 



and 



Fo(x) 



EVC^M/mjj) 



1 — Qx + 3x 2 + 2x 3 — 6x 2 \nx 



(13) 



(14) 



(15) 
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Since M& < m should be satisfied for N^ dark matter, the function F 2 (xk) can vary only 
between 1/12 (a:* = 1) and l/6(x^ = 0). To suppress the branching fraction B(fi —>■ &y) 
which is inversely proportional to the fourth power of mo, we need a large value of mo. 
On the other hand, the observed dark matter relic density Vtji 2 ~ 0.12 requires m to be 
below 350 GeV for yk = (Y.a/3 \hakh*n k \ 2 ) 1 ^ ^ 1. This means that if | Ya h^h^] appearing in 
Eq. ((H]) is also of order 1, then B([x — > ej) ^5x 10~ 7 , which is several orders of magnitude 
above the experimental upper bound of 1.2 x 10~ n . 

To satisfy the /i — > cy constraint, we consider the possibility that M\ t 2,3 are nearly 

degenerate. In that limit, 

(m x e~ iai 0\ 



{M v U= I{M2 ^ ) Y.hc A h pi = U* 



m 2 e- ta2 

V m 3 J 



U\ 



(16) 



where U is given by Eq. (JSj). A simple solution is then 



Mm; 



I(M 2 /ml) 



1/2 



e -ta,/2 fj* with = 



(17) 



Then we obtain 



C 2 



F 2 (M 2 /m 2 )M 
I(M 2 /ml) 

F 2 (M 2 /m 2 )M 
I(M 2 /ml) 



i 
S12C12 



(m 2 - mi) + 



size 



-{c\ 2 m 1 + s\ 2 m 2 -m 3 ) 



Thus the suppression of C 2 is possible because m 2 — mi is related to Am^, and c\ 2 mi + 
s\ 2 m 2 — m 3 is related to Am 2 ltm in neutrino oscillations jTHllIlj. 
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Figure 4: C 2 versus mi in the case of normal ordering for S13 = 0.1 (solid), 0.05 (dash), 0.01 (dot-dash), 
where C 2 = 4.6 x 10~ 4 (horizontal line) corresponds to the experimental upper bound B(/i — > e~f) — 
1.2 x 10- u . 



Let us assume 5 = and consider the normal ordering of neutrino masses, i.e. m 3 
is the largest mass. We then set h 3 = (J2 a \h a 3\ 2 ) = 1 which is equivalent to having 



Mm 3 /I(M 2 /ml) = 1. Hence 



C 2 



0.067\ 



m 3 



J |ci 2 (si 2 - s 13 ci 2 )(m 2 - mi) - s 13 (m 3 - m 2 )\ < 4.6 x 10" 



(19) 



where F2 = 0.0948 (corresponding to mo = 345 GeV and M = 290 GeV). Using Am^ = 
Am 2 sol = 7.9 x 10 -5 eV 2 and Am§ 2 = Am^ m = 2.3 x 10 -3 eV 2 , we plot C 2 versus mi 
in Figure 0] for S13 = 0.1,0.05,0.01. The horizontal line is the experimental bound C 2 = 
4.6 x 10~ 4 corresponding to B(/j — > eq) = 1.2 x 10 -11 . We find that for s n k, 0.26, this 
constraint cannot be satisfied. For S13 less than its experimental bound of 0.2, there is a 
lower bound on mi according to the approximate empirical formula 



mi 
eV 



^ 0.02 + 1.4|si 3 - 0.02| - 2.9|si 3 - 0.02| 



(20) 



except for a tiny region near mi = and S13 = 0.09, and a small region near mi = 0.01 eV 
and S13 = 0.04. In Figured we can see that the plot for S13 = 0.1 is getting close to the first 
region from its dip at vrt\ = 0, and that the plot for S13 = 0.05 has a small allowed range 
near mi = 0.01 eV. 
O.l 
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Figure 5: C 2 versus 7713 in the case of inverted ordering for S13 = 0.2 (solid), —0.05 (dash), 0.0 (dot-dash). 



In the case of inverted ordering, i.e. m 2 is the largest mass, we are already guaranteed 



that mi : 2 > JAT7132 ~ 0.048 eV. For completeness, we set hi = 1 and plot C 2 versus 777,3 i n 
Figure for s 13 = 0.2,-0.05,0.0. Here, for s 13 <: 0.24 and Si 3 ^ —0.27, the experimental 
constraint cannot be satisfied. In other words, the constraint on #13 from \x — > cy coincides 
roughly also with that from neutrino oscillations. 

For the simple solution of Eq. (16), as we can see from Figure 4 and 5, all the neutrino 
masses may be assumed to be degenerate to satisfy the /1 — > e-y constraint; hence the effective 
mass (m ee ) in neutrinoless double beta decay is approximately given by 

(m ee ) ~ 777x10.572 + 0.428 cos(«i - a 2 )| 1/2 . (21) 

We also allow the heavy iVj masses M lj2j3 to be slightly different, so that our approximation 
that only one of them is the candidate for dark matter remains valid. Note that AM/M 
only needs to be of order 10~ 3 for Eq. (18) to be valid. In analogy to [i — > ej, there are 
also dark- matter contributions to r — > fj,j, r — > ej, and the anomalous magnetic moment 
of the muon. However, they are at least one or more orders of magnitude below the present 
experimental bounds. 

In conclusion, we have shown how cold dark matter and neutrinoless double beta decay 
may be connected if neutrinos acquire mass only because of their interactions with dark 
matter. We repeat the basic argument presented earlier. The existence of dark matter 
requires a class of new particles which are odd with respect to an exactly conserved Z 2 
symmetry. Their interactions with neutrinos and charged leptons must not be too weak to 
be identified as cold dark matter with the correct value of their relic density in the Universe 
at present. On the other hand, they are also responsible for the masses of neutrinos and their 
observed mixing in neutrino oscillations. This implies necessarily flavor changing transitions 
such as fj, — > ej. In order to suppress the latter, the parameter space of neutrino masses is 
limited, thereby enforcing a lower bound on neutrinoless double beta decay. For N as dark 
matter, this is typically of order 0.05 eV, even though much lower values are still allowed from 
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accidental cancellations. More importantly, this connection between cold dark matter and 
neutrinoless double beta decay can be tested in the near future at the Large Hadron Collider 
and complemented by a host of experiments on neutrino oscillations and neutrinoless double 
beta decay already under way and being planned. 
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